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Abstract 

A path integral is presented that solves a general class of linear second order par- 
tial differential equations with Dirichlet /Neumann boundary conditions. Elementary 
kernels are constructed for both Dirichlet and Neumann boundary conditions. The 
general solution can be specialized to solve elliptic, parabolic, and hyperbolic par- 
tial differential equations with boundary conditions. This extends the well-known 
path integral solution of the Schrodinger /diffusion equation in unbounded space. 
The construction is based on a framework for functional integration introduced by 
Cartier/DeWitt-Morette. 
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1 Introduction 



Functional integral solutions of second order elliptic and parabolic partial differential 
equations (PDEs) have been known for a long time. In particular, stochastic methods 
have been used to solve inhomogeneous elliptic and diffusion-type parabolic PDEs 
(see e.g. [1], [2], [3]), and path integrals ""^ in physics have been used to solve the 
Schrodinger equation ([4], [5]) and homogeneous Dirichlet problem 



Stochastic methods have perhaps been the most fruitful in terms of solving PDEs. The 
majority of stochastic results have been obtained for scalar functions on M"; although 
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Elworthy ([7]) was able to extend some of these to tensor fields on manifolds. However, 
stochastic methods are limited in scope because they are based on the Wiener process. 
Consequently, they cannot be applied to Schrodinger-type parabolic PDEs (except 
through analytic continuation), and they have nothing to say about hyperbolic PDEs. 

On the other hand, path integrals have been developed extensively in a myriad of 
physics applications. They have been extended to tensors on manifolds (following 
Elworthy), Grassman variables, and fields. However, for the most part, path inte- 
grals have not been exploited for solving PDEs outside of Schrodinger-type parabolic 
equations in unbounded space. 

In this paper, a general path integral is presented that solves a wide class of linear 
second order PDEs with given Dirchlet/ Neumann boundary conditions. Relevant ele- 
mentary kernels are also constructed to facilitate incorporation of Dirchlet/Neumann 
boundary conditions. In a subsequent paper, the general path integral will be special- 
ized to solve elliptic, parabolic, and hyperbolic PDEs; several examples will be worked 
out to check the validity of the constructed solutions against known solutions; and 
some new calculational techniques will be introduced. 

Section 2 contains the outline of a framework for functional integration developed 
by Cartier/DeWitt-Morette in [5] (see also [8]). It allows one to define path integrals 
in a general setting. The three main ingredients of the framework are: an infinite 
dimensional space of pointed paths which take their values in a (complex) manifold; 
a parametrization of this infinite dimensional space of pointed paths by a Banach 
space; and integrators (that assume the role of problematic measures in infinite di- 
mensional spaces) defined on the parametrizing Banach space. The scheme is then 
amended by restricting the space of integrable functionals; and a new integrator, the 
gamma integrator, is introduced. These additions play a key role in the construction 
of the solution of the general PDE. Speciahzing various components of the functional 
integration framework leads to path integrals that are solutions of the three classes 
of PDEs mentioned above. 

Part of the utility of the Cartier/DeWitt-Morette scheme is that it deals with inte- 
grators instead of measures. For example, the integrator corresponding to the Weiner 
measure can be used to solve diffusion-type parabolic PDEs. However, while it is 
not possible to define a measure that corresponds to the complex version of a Weiner 
measure, it is possible to define an integrator that does. Consequently, in this scheme, 
one can solve Schrodinger-type parabolic PDEs directly. 

Since the Cartier/DeWitt-Morette formulation includes parametrizations that are 
similar to stochastic differential equations and integrators that are similar to Weiner 
measures, it is not surprising that some of the solutions constructed here are similar to 
the stochastic constructions. Some of the similarity arises from the use of a dependent 
variable that reparametrizes the time along a path. From a physics perspective, the 



2 



variable is a non-dynamical degree of freedom that is introduced because the paths of 
interest are point-to-boundary instead of point-to-point — the presence of the bound- 
ary induces a constraint. Of particular importance is the minimum time to reach a 
given boundary starting from some point. Alternatively, one can consider the path- 
dependent time to reach a given boundary starting from some point. In stochastic 
parlance this variable is the first exit time from the boundary (the first time at which 
a given stochastic path intersects the boundary). Because these two notions are simi- 
lar, they lead to similar constructions. However, there is a crucial difference that will 
become clear later, and the construction based on the former notion is advantageous 
in terms of evaluating path integrals. 

Section 3 constructs the path integral solution along with its associated kernels for 
various boundary conditions. The gamma and Gaussian integrators play a major 
role in the construction. The reader should be warned that no attempt is made here 
to determine conditions for existence and uniqueness of the constructed solutions. 
However, it is reasonable to expect that the path integral setting will be a valuable 
tool in addressing these difficult issues, because it shifts the focus to the function 
space instead of the target manifold. 

One of the virtues of path integrals is that they automatically incorporate boundary 
conditions; sometimes implicitly, sometimes explicitly. As such, they offer an alternate 
approach to study and solve PDEs with non-trivial boundary conditions. In particu- 
lar, complicated geometries are no different from symmetrical geometries in principle. 
Since computational and approximation methods of path integrals are fairly well de- 
veloped, finding solutions for a given boundary problem should be at least systematic 
if not simplified. 

Supporting details of the functional integration framework are presented in Appendix 
A and Appendix B. Appendix A contains several properties of functional integrals 
that are fairly straightforward but many of which have not appeared explicitly in the 
references. Appendix B discusses the well-known Gaussian integrator and introduces 
the new gamma and Hermite integrators. The Hermite integrator is a generalization 
of the Gaussian integrator and the gamma integrator is the scale invariant analogue 
of the (translation invariant) Gaussian integrator. 



2 Path Integration 

2.1 Cartier /DeWitt-Morette scheme 

According to the general scheme ([5], [8]), a path integral is defined on a separable 
Banach space X with a norm ||x|| where x e X is a map x : E — > M. E is a 1- 
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dimensional manifold and M is an m-dimensional manifold. The dual Banach space 
X' 3 x' is a space of linear forms such that {x', x) e C with an induced norm given 

by 

= sup I {x', x) |/||a;|| . 

Assuming X' is separable renders it a Polish space, and, therefore, it admits complex 
Borel measures /x. 

Given two continuous, bounded, and /x-integrable functionals Q : X x X' ^ and 
Z : X' ^ C''', an integrator on X is defined by 



Definition 2.1 



/ e{x,x')Ve,zx:^Z{x') . (2.1) 



A space, ^{X), of functionals integrable with respect to this integrator consists of 
functionals defined by 



Definition 2.2 



F^(x) := / e{x,x')dii(x') . (2.2) 



In this paper, F^{x) for some path x E X will be a tensor field or tensor distribution of 
type (r, s) along xiT,) on the manifold M. The tensor distributions (suitably restricted) 
will be used to construct kernels of second order linear partial differential operators. 
In the sequel, all tensor fields and distributions are assumed to be well defined, i.e. 
they are continuous, differentiable, integrable, bounded, etc. as required for any given 
case. Since existence and uniqueness issues are not addressed, this simplification is 
warranted. 

The map /x 1-^ is required to be injective rendering T{X) a Banach space, and 
J-{X) can be endowed with a norm defined to be the total variation of ji. An 
integral operator T^e,zx the normed Banach space J-{X) is then defined by 



Definition 2.3 



[ F^ix)Ve,zx:= [ Z{x') dfi{x') . (2.3) 



The integral operator zx is a bounded linear form on J-{X) with bound |/^ Vq zx\ < 

\\F,\\ ([8]). 

A fruitful generalization is to consider functionals and Z that depend on a param- 
eter A where, for example, A may be in R or C. It will be assumed that F^(x; A) is 
differentiable with respect to A, and that 9^"^F^(x; A) e J^iX) for n e {1, 2}. 

The right-hand side of (2.3) should not be thought of as a prescription for calculat- 
ing the left-hand side — in practice, one neither knows nor specifies /i. Instead, one 
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invariably makes use of some form of localization to reduce the left-hand side to a 
finite dimensional integral. Hence, a 'good' (or at least useful) characterization en- 
coded in Definition 2.1 must reduce to the correct finite integral for any dimension: 
consequently, it is reasonable to expect path integrals to possess properties analogous 
to finite dimensional integrals. Some properties of path integrals that are particularly 
relevant are listed in Appendix A. 

For actual applications, the goal is to define an integral over an infinite dimensional 
space, -PJM, composed of pointed L^'^ maps x : E — > M with a mutual fixed end-point 
Xa e M. Here, M is a (possibly complex) m-dimensional paracompact differentiable 
manifold, and I will restrict to S = I C R or E = 5'"'" and refer to x E PjM as a 
path (hence the term 'path' integral). Additionally, assume a set of c? < m linearly 
independent vector fields X(q,) where a e {!,..., d}. The set {X(q,)} generates a 
sub-bundle V C TM. 

However, in general PjM. is not a Banach space. Analogous to the definition of inte- 
gration on finite dimensional manifolds in terms of integration on M", path integrals 
on PJM can be defined in terms of path integrals on a separable Banach space. The 
definition relies on a parametrization of the paths effected by the differential system 

dx(t,z)-Y{x{t,z))dt^:^^^){x{t,z))dz\t) 

(2.4j 

2^(to) ~ 

where Z 3 z is & Banach space, and a;(t) — Y(x{t,z)) e Yx{t) Q T!r(t)M for each 
t e E and some given vector field Y on M. The solution of (2.4) will be denoted by 
a;(t, z) = Xa ■ S{t, z) where Z'(t, z) : M ^ M is a global transformation on M such 

that X{ia, ■) =^a- ^(ta, ") = X„. ^ 

The nature of the integrator on Z dictates the range of the paths z e Z (or vise- versa). 
For example, for a translation invariant integrator, z G PqW^{'C^) =: Z^; and for a 
scale invariant integrator, z G P^]R+(C+) = In(Z^) =: Z+ with z(ia) = for some 
t(j e E.^'^'^ The differential system (2.4) associates a pointed path z G Z^{Z^) 
with each pointed path x G PjM yielding a parametrization P : Z^{Z^) — > PjM. 

The parametrization P : Z^{Z^) PjM and Definition 2.3 allow for a rigorous 
definition of the path integral of a functional F{x; t) over PjM: 



^ Compactness of the manifold M may be required in some cases to ensure global existence 
and uniqueness of the map z ^ x. However, one is usually interested in a submanifold 
U C M so compactness of the entire manifold may not be necessary. 
3 The notation P^R'^{C'^) is short-hand for P^R'^ or P^C^. 
^ C+ denotes the right complex plane. 

^ Translation and scale invariant integrators are characterized in Proposition A. 5. 
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Definition 2.4 



^.^ Restrictions 

In order to implement the analogues of invariant measures and integration by parts 
in finite dimensions, it is useful to impose two restrictions on the general scheme 
outlined in the previous subsection. 

Let y be a separable Banach space and M : X — F be a diffeomorphism with 
derivative mapping M^'^-, : T^X — TyY. The derivative map possesses a transpose 

W' -.TY' ^ TX' defined by 

(M\y,)y', x)t^x := {y', M[,^x)T,y . (2.6) 

M induces a germ R .Y' ^ X' along with its associated derivative mapping -R'^^^/) : 
TyiY' — > Tx'X' and transpose R\x) '■ TxX — > TyY from the relation 

{R[^,^y',x)T^x^{y',RUx)T,Y. (2.7) 

Evidently for y = M{x) and x' = R{y'), M' = R' and R! = M' . 

Define the /x-integrable functional Q -.Y xY' ^ C('^+") and 1 -.Y' ^ C('^+") by 

eo (M X i?-^) := e (2.8) 

and 

Zoii:-^:=Z (2.9) 
along with their associated integrator 

/ Q{y,y')V^;gy:^Z{y'). (2.10) 

Let R{Y') = Y'. The first restriction relates Z and Z at the same point in Y'; 
Restriction 2.1 

Z{y')^\DetR[y,)\-'Z{y') (2.11) 

for R a diffeomorphism and non- vanishing |Deti?|j^,-) |. For the determinant to be well 
defined, require that R' be nuclear (see for example [15] for relevant details). 
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To state the second restriction, it is useful to interpret Ffj,{x)T>Q^zx =: Ffj_ as a form 
on X. Let J-'^{X) denote the space of integrable forms on X. Let dy denote the 
Gateaux derivative of Fi^{x) in the y direction, i.e.. 

Assume F^{x) and Vq^zx are Gateaux differentiable, then require 
Restriction 2.2 



/ dyF^ = [ dy[F^{x) Ve,zx] = (2.13) 



for arbitrary y ^ X. 



These two restrictions have seemingly been "pulled out of a hat", but, as shown in 
Appendix A, they enable the characterization of invariant integrators and integration 
by parts. Henceforth, JF^(X) will denote the space of integrable functionals restricted 
by (2.11) and (2.13). Likewise, ^r{X) is assumed to be similarly restricted. Note that 
(2.11) and (2.13) may be quite restrictive and J^r{X) may be severely hmited or even 
empty. ^ 



3 Path Integral Solution of PDEs 



The path integral framework introduced in Section 2 allows the construction of solu- 
tions to general classes of linear second order partial differential equations. It turns 
out (as shown in a subsequent paper) that various choices of the manifold M and the 
parametrization P will yield elliptic, parabolic, and hyperbolic PDEs. It is noteworthy 
that the construction can cover all three cases. 



3.1 General Solution 



It is now possible to construct a path integral for PDEs with finite boundary. The con- 
struction depends crucially on a time reparametrization r and its associated integrator 
Vt. Appendix C presents the motivation for introducing a time reparametrization for 
bounded regions. It should be pointed out that the reparametrization used here is fun- 
damentally different from the time reparametrization sometimes employed in stochas- 
tic methods (sec e.g. [17]) or by Duru and Kleinert ([18]): there the reparametrization 
is path dependent. Specifically, stochastic constructs parametrize the first exit time of 

® Clearly, the question of existence of PDE solutions constructed in the next section is 
directly related to the nature of J^r{X). 
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a given path — since the paths depend on a random variable, then so docs the first exit 
time. Here the reparametrization is path independent and the first exit time used in 
this construction is for a critical path. This difference allows z and r to be integrated 
independently. 

With these issues in mind and some refiection on existing functional and path in- 
tegral solutions to selected PDEs, a path integral solution of a linear second order 
inhomogeneous PDE is constructed. 

First some definitions: 

Definition 3.1 T+ is the space of L'^'^ functions t -.T — [ta, — >■ [to, Tb] U [ra, r^]* C 
C+ such that \dT/dt\ > and Ta '■— T{ta) — 0. 

PJM is the space of L^'^ pointed paths x : r(T) M with x{Ta) = x^ G U C M whose 
velocity vectors are elements o/ V C TM. 

is the space of L"^'^ pointed paths z : t(T) R'^{C'^) with z(t„) = 0. 

Definition 3.2 Letx{T{t),z) = Xa- S{r{t), z)— where i:{r{t), z) : M ^ M is a global 
transformation on M such that a;(r(ta), z) — x^ - S{T{ta), z) — x^ — denote the solution 
of the parametrization P : ^ PjM according to 

' dx{T{t), z) - Y{x{T{t), z))dT = X(,)(x(r(t), z))dz''{T{t)) 



Definition 3.3 The functional S{x{Tb, z)) is defined by 

S{x{n, z)) := 7rQ{x{n, z)) - T Vix{t, z)) dt (3.2) 

Jo 

where 



rn rn 

Q(x(Tb, z)) := (Dx, x) = / / xHt, z)Dij(t,'u.)xHu, z) dtdu (3.3) 

Jo Jo 

in a local chart. 

The operator matrix Dij{t,u) with i,j G {!,..., n} depends on the nature of M. 
In physical applications, Dij is usually related to a metric or symplectic form. For 
example, Q is often identified with an action functional of the form 

A{x) = r /i,(t) (i(t) - Y{x{t)), x{t) - Y{x{t))) dt (3.4) 

where hy_ is a quadratic form on TxM and Y e TM is some vector field. The end- 
point value of the Green's function of D is the matrix G"^ :— G°'^{tb, %) discussed in 
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Appendix B.l. 



Definition 3.4 Let denote the first exit time of a critical path, with respect to 
S{x{ti), z)), starting at a pointXa G U C M; i.e., Xcr{T^J ='■ e 9U and Xcr(O) — x^. 

Recall from Appendix C that a critical path for a variational problem with a boundary 
must satisfy the Euler equations and supplemental 'transversality' conditions. 

Definition 3.5 Let Q := X T+. Characterize an integrator on Vl by := 
Vq wzVt where Vq ^z is a Gaussian integrator characterized by (B.2) and Vt is 
the gamma integrator characterized by (B.55). 

Theorem 3.1 Let Ml be a real(complex) m-dimensional (m > 2) paracompact dif- 
ferentiable manifold with a linear connection, and let V be a bounded orientable open 
region in M with boundary dV. ^ Let f and ip be elements of the space of sections 
or section distributions of the (r, s)-tensor bundle over M. ^ Assume given the func- 
tional S{x{Ta' , z)) whose associated bilinear form Q satisfies Re((5(a;(ra', 2;)) > for 

{x{Ta>,z))^Q. 

Ifxi^a ■ S{{t',t),z)) e J^nin) where' 

Xi^a- ^{{r',r),z)):= / 0{{r',r) - t„/)/(x„ • i:(Ta/, 2;)) exp {-5(x(Ta/, 2;))} dra' 

+ ^{xa - S{{r',r),z))exp{-S{x{{r',r),z))} ; (3.5) 
then, for Xa = x{Ta) G U, 

S^(x„) = x(x„ • r(ri, z)) VQ (3.6) 
is a solution of the inhomogeneous PDE 



with boundary condition 



4^(x„) = -/(x„) (3.7) 



^{xb) = ip{xB) . (3.8) 



^ The issue of regular and irregular points of 5U will be ignored. However, I assume sufficient 
regularity of 9U when required. 

^ All tensors and tensor distributions are assumed to be well defined, i.e., they are contin- 
uous, differentiable, analytic, bounded, integrable, etc. as required for any given case. 
^ From a physics perspective, it may be useful to express % as a functional of r and add 
a term to S{x{t, z)) that restricts r to the value r^. However, for our purposes, this is an 
unnecessary step. 
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The symbol >Cx(„) represents the Lie derivative in the X(a) direction. (Note that it 
is possible to have highly non-trivial functions of x pre-multiplying partial derivative 
terms due to the presence of Lie derivatives.) My convention for the step function of a 
complex variable is ^^(|t| — |to|) = if |t| — |to| < and 6'(|t| — |to|) = 1 for |t| — |to| > 0. 
A specific choice of integration path C+ C C+ is dictated by a particular application, 
which usually entails a restriction on r G T+. For example, in typical applications 
either r = r* or r = — r* so that C+ C or C+ C iR respectively. 



Proof of Theorem 3.1. From the definition of 3^(xa) and xi^a • '^('^i? ^)) Propo- 
sition A. 3, 



Lfl^(x„) :-- 



where 



Qa/3 

An 



(3.9) 



{UrJ){^a) ■■= / /(x„-r(v,^))exp{-cS(x(v,z))} Vq^wZ. (3.10) 



By Proposition A. 5, the integrals over T+ can be reduced to one dimensional integrals 
over C+. The specific form is given in (B.56). Then assuming the necessary conditions 
on the integrands, L can be taken inside the integrals. 

Furthermore, 

Lemma 3.1 L commutes with U^^t-',t}- 
Proof. Consider 

{U^r',r}^){^a) = f ^ Cp{x{{t' , t) , z) BXp {-S{x{{r\ t) , z))} Vq,wZ . (3.11) 

Taylor expand ifi{x{{r', r), z) and cxp {—S{x{{r', r), z))} (assuming the necessary dif- 
ferentiability) about X(i, and use Proposition A.l to interchange the sum and the inte- 
gral. Take the terms that depend explicitly on x^ outside the integral over Z^, apply 
L, and reverse the steps. (Recall that the resulting functional have been assumed 



10 



integrable in Section 2.) □ 



Lemma 3.2 



/ e{{T\ t) - Ta')L{UrJ){Xa) dTa' = -/(x„) - (?7(.,,.)/) (x„) . (3.12) 



Proof. Scale the time variable by {t',t)~^ in the parametrization (3.1). The space 
of paths gets mapped to the space of paths with elements z(r(t)/(r', r)) = 
(r', T)~^/^2;(T(t)), and the parametrization becomes 



i(r(t)/(r',r))-(r',r)Y(x(r(t)/(r',r))) = 

(r', r) V2X(„)(a;(r(t)/(r', r)))2 "(r(t)/(r', r)) 

x{Ta/{T',T)) =X„ . 



(3.13) 



Expand f{x{{T',T),z)) about the point Xa making use of (3.13): 



/(x„ • U{{r',r),z)) = f{xa) + (r', t)V2£x(„/(x„)?°(1) 

+|(r',r)£x(„/x(,,/(x„)?'^(l)?^(l) 
+ (r',r)£Y/(x.) + 0((r',r)3/2). 



(3.14) 



Likewise, the functional S{x{{t',t),z)) becomes 



S{x{{t',t),z))=7tQ{x{1,z))-{t',t) V{x{t,z))dt (3.15) 

Jo 



and 



y(x„ • S{{r',r),z)) = y(x„) + (r', r) V2£^^^^y(x„)?°(l) 

+1 (r',r)£x(„,£x(,,y(x„)z"(l)z^(l) (3.16) 
+ (r',r)£Ynx„) + 0((r',r)3/2). 



From Appendix B.l — which records the normalization, mean, and covariance of a 
Gaussian integrator — and (3.14)-(3.16), 
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+l(r',r)£x(„)>Cx(,)/(x„)^"(l)?'^(l) 
+ (r',T)£Y/(x„) + 0((r',r)=^/^)} 
x{l + (r',r)\/(xJ + 0((r',r)2)}l).;(J) 
= /(x„) + (T',T){L/}(x„) + 0((r',r)3/2) . (3.17) 

The last equality is a consequence of Proposition A.l. It follows from Lemma B.2 
that {U(^r',T)f){^a) satisfies the PDE 

with initial condition 

(^r',r)=0/)(Xa)=/(Xa). (3.19) 

The statement of the lemma is the integral form of the PDE. □ 
Using Lemmas 3.1 and 3.2, equation (3.9) can be written 

Lfl^(x<,) = -/(x<,)-A^/ (C/t{/-Lv?})(x„)d(lnt) . (3.20) 
0= / (^7t{/-L<^})(x„)d(lnt) . (3.21) 



Lemma 3.3 



Proof. By Proposition A. 3, 



= 4 X, X(xa ■ Sirt, z)) Vq,wzVt 

(f7.xX)(x„)PT. (3.22) 



From Proposition B.9, Vt is scale invariant and it therefore verifies the integration 
by parts formula (A. 20) for t/^x x ^ ^r{T^)i 

= / ^ °^ Vr . (3.23) 

But, 
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ST{t) 



-ill K{r\ r) - Ta')f{x, ■ S{Ta', z)) exp^'^ ^^^^<^"^m dr,, 
Jzi \Jc+ 

+ b'(x„ • E{{t\ t), z)) - <^(x„ • E{{t\ t), z))S'{x{{r', r),z))] 

X exp^-'^^^^^"''"^'"))^} Vq^wz 

= / {-/(x„.r((r',r),z)) 

+ [<^'(x„ • S{{r',T),z)) - <^(x„ ■ E{{T',T),z))S'{x{{T',r),zm 

X exp{-'5(^(<^''^>'^»> Vq^wz 

/ {-/(x„ • S{{T',T),z))+^{xa ■ i:((r',r),^))}exp{-^W<-'^)'^)» Vq,wz 

= -(^(r',r>{/-$})(Xa). (3.24) 



From Proposition B.2, 

U{r',T){f - ^} = t^(T',T)+(5t.,T){/ - ^} = t^(r',r>(t^(5t.,T){/ " ^}) 



(3.25) 



since {St^^r) = 0. Use the expansion technique from Lemma 3.2 to expand (p' and S' 
in do the integration over as before, and evaluate at {St^ , r) = to get 



W-^}(Xa)=/(Xa)- 



^v + ^/:x(„,/:x,,, + nx) 



V(Xa) 



:/(x„) - L(fi{Xa) . 



(3.26) 



Hence, (3.23) can be rewritten as 



= ^^(C/,x{/-L<^})(x„) Vr 



(3.27) 



Since this reduces to a one dimensional integral according to (B.56), the lemma is 
verified. ^° □ 



So, by Lemma 3.3, equation (3.20) reduces to the inhomogeneous PDE 



(3.28) 



"'^'^ Intuitively, this makes sense — at least in the case where {Ut{f — L(p})(xa) solves the 
diffusion equation, i.e. t G M_|_. The integrator is scale invariant, and {Ui{f — L<^})(xa) — > 
as t — >■ cxD. 
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To verify the boundary condition, note that the definition of imphes that x, 
Xfi e dV imphes 0. The boundary condition foUows readily: 



-a 



fl^(Xa)| 



Xa=XB — 




v(xb) . 



(3.29) 



The second equahty holds because ■ S{0,z) = xb is independent of z, and the 
third equality follows from the normalization of the integrators. 

This completes the proof of the theorem. □ 
Remarks about Theorem 3.1: 

— Note that there may be more than one critical path and, hence, more than one r^^. 
In that case, (3.6) should then include a sum over all r^. 

— If the boundary of U is at infinity and ^(xb) ~ 0, then (3.6) can be written as 
the Fourier/ Laplace transform of a path integral which solves an inhomogeneous 
parabolic equation associated with (3.7). This is a well known result for the case 
of elliptic PDEs. 

However, when M corresponds to the phase space of a dynamical system, the path 
integral is more naturally rewritten in terms of a Lagrange multiplier by employing 
the map i? : r+ — > yl which maps r i— > A by 



Consequently, the Fourier transform interpretation is no longer evident. 

Under this map, the action then takes the form of a constrained system and the 
domain of integration reverts back to [t^, t;,]. This justifies expressing (3.6) in terms 

of a path integral over Tj~ instead of an equivalent but simpler integral over C+. In 
this sense, the phase space representation of (3.6) appears to be more general and, 
therefore, perhaps more useful. 
- For f — 1, (f = 0, and suitable V{x), 



This is easily seen, for example, by using normal coordinates near the boundary and 
recalling that critical path. As soon as the starting point lies on the boundary, 

= for any endpoint on the boundary. 

If the potential F(x) has no non-vanishing constant term, then the associated parabolic 
PDE has a different potential, viz. F'(x) = V"(x) — 2mS. 




(3.30) 
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so that L{t^Jj,+ — —1 with boundary condition (r^)^+ = 0. (More precisely, this 
follows after integrating (3.31) over Xa'-) 

— The (pointed) path integrals can be transformed into (pointed) loop integrals by 
replacing the parameter interval T = [ta, t;,] with and requiring r to be a home- 
omorphism. 

— The case of simply connected, compact U without boundary can be included in 
the theorem by taking M to be a simply connected, compact manifold without 
boundary and then deforming the boundary dV to a point (without intersecting 
Xa, the point of interest). This amounts to putting cp = 0, and letting r,^ represent 
the period of a critical loop in U. 

— It is important to note that the proof of the theorem doesn't require that dV 
be closed: Nor does it address the boundary type and boundary conditions nec- 
essary for existence and uniqueness of the solutions. For scalar functions on M", 
typical parabolic and hyperbolic PDEs require an open boundary with Cauchy 
initial condition — possibly supplemented with Dirichlet/ Neumann boundary con- 
ditions. On the other hand, eUiptic PDEs require a closed boundary with Dirich- 
let/Neumann boundary conditions. However, the boundary type and boundary 
conditions necessary for existence and uniqueness of solutions for the generalized 
PDEs considered here are open questions. 

3.2 Kernels and Eigenf unctions 

In practice, finding the functional x for general boundary conditions can be difficult. 
The task can be simplified by constructing kernels of (3.7), denoted by K, appropriate 
for Dirichlet/Neumann boundary conditions. 

3.2.1 Dirichlet Kernels 

Lemma 3.4 The Dirichlet elementary kernel is given by 




xexp{-S{x{Ta',z))} dTa'Vn (3.31) 

where Xa,Xa' G U and 6{xa,Xa') denotes a Dirac bitensor composed of Kronecker delta 
symbols ( which are collectively denoted by 1 since the specific form is not necessary 
for our purposes) and the scalar Dirac delta function 5(xa,Xa/). 

Proof. It follows immediately from the theorem that -R'[f is the elementary solu- 
tion of the inhomogeneous PDE. Moreover, K^f ^ (xB,Xa/) = since x^ — > x^ implies 
— > 0. Evidently, ■'(xa,x<j/) is the elementary kernel of the PDE which vanishes 
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on the boundary. □ 



Lemma 3.5 The Dirichlet boundary kernel is given by 

Kf\x,, xs) := 5{x{t^,z),xb) exp [-S{x{t^^, z))} VQ (3.32) 
where xb G dV. 

Proof. According to the theorem, K^^\xa,XB) is a kernel of the homogeneous PDE 
and K^^\xb',xb) — 5(xb/,xb) for x^' G dU. □ 

Corollciry 3.1 For Dirichlet boundary conditions (^{xb), the solution of the inhomo- 
geneous PDE can be written 

«^(X„) = / iir{f)(x„,X„0/(x„0 dXa' + I Kf\xa,XBMxB) dxB (3.33) 

where dxa is the volume element on V, and dxB is the associated Leray form on dV. 
Proof. The corollary follows from the preceding lemmas. □ 



Notice that the boundary kernel is a path integral over paths starting at x<j at and 
ending at some point on the boundary at time t^^. It can be expressed indirectly in 
terms of paths starting at x^ at Tq and ending at the specific point Xcr{T^J =: x^ e dU. 

This follows because K^^\xb',xb) — S{xb',xb) for xb' G dU. (The point is it doesn't 
cost anything to propagate along the boundary.) Hence, 



) (X„, X^) = J^^ Kf^ (X„ Xb)K^^^ (xb, X^) dXB 

= / Kf\xa, XB)d{xB,X^) dxB 

dixiTi,z),xi) exp {-SixiTi,z))} Vn 



(3.34) 

This form is particularly useful in calculations when one expands about the critical 
path(s). 

For later use, it is convenient to rewrite (3.31) for U embedded in a non-compact M 
for the special cases when r = r* or r = — r*. 
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Proposition 3.1 

lif if ^ (x„ ,^a>)^ LS{X {Ta' , z) , X„/ ) exp { (t„/ ,z))} dTa> Vfl 

JCl UM. 

- / LS{x{Ta',z),Xa')exp{-S{x{Ta',z);s)}dra' Vn 
Ja Jr 



(3.35) 



where x{Ta' , •) = a{xa)-S{Ta', •) ; cr(xa) = Xcr(T^); anc?M (denotes R+ oriR depending 
on whether t — t* or t — — r*. 

Proof. For r = r*, 



K'^\xa,Xa') = f I 6{x{Ta>,z),Xa>)eX^{-S{x{Ta>,z))} dTa' 

Jn Uo 

^ <5(x(ra', z),Xa/) exp {-5(a;(ra', 2;))}(ira 
-K'i(xa,Xa') - F^(x„,x„/) ; 



(3.36) 



and for r 



-R^lf^(Xa,Xa') = [iiL+ (x„,Xa/) - FJ(x„,X„/)] + [K^ (x„, X„/) - Fu(Xa,X„/)] , 



(3.37) 



where 



-K"oo(Xa,Xa') -Fu(x„,X„/) 
- rO 

/ S{x{Ta', z),Xa>) exp {S {x {Ta' , z))} dTa' 
.J —oo 



j 

Jn 



/ " d{x{Ta' , z),Xa>) exp {-S{x{Ta' , z))}dTa' 
J —oo 



(3.38) 



Let iifoo denote K'^ or [JT^ + K'^] and Fu denote F^ or [Fj + F^] depending on 
whether Imr = or Rer = respectively. According to the theorem, K^oi^ay^a') 
is the elementary kernel for the boundary at infinity since x^ — oo implies |Tj^| — >■ 
oo. ^'^ It follows that Fu must be the kernel to the homogeneous PDE with boundary 
condition 

Fu{xB,Xa') = K^{xB,Xa') . 



More precisely, since the paths are L^'^, their energy, and hence velocity, are bounded. 
This requires Ix,;^ | — > oo as Xa — > oo. 
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(Roughly speaking, Fu satisfies the homogeneous PDE because Xcr{Ta') ^ U for Ta' > 
<■) 

Fu can be rewritten by introducing a point transformation cr : U — > U by i— > 
cr(xa) = XcriT^J — (which is a point on dU). Consequently, r,^ i— > t^^^^^ = under 
this transformation. 

In order to leave the boundary conditions unchanged for paths whose end-points are 
restricted by the delta functional in Fu, it is necessary to simultaneously transform 
Xa' Xa' such that X(j/ = a~^{xa'). To see this, let S denote the transformation from 
Xa to the final point x^/ and write 

•) = a{Xa) ■ S{Ta', ■) =Xa- S{Ta', • ) = ^a' (3.39) 

SO that S = S o a. Then, for paths with transformed initial point a{xa) and fixed 
end-point x^', 

cr{x{Ta', ■ )) = a{xa') = Xa' = x(v, • ) = x„ • ^(v, • ) . (3.40) 

In particular, 

^{xUrt, ■))=^i (3.41) 
remains consistent with the the parametrization and boundary conditions. 

For paths that contribute to Fu, it follows that 

Fu{a{xa),a~'^{xa')) 

= / LS{x{Ta',z),Xa')exp{-S{x{Ta',z);s)}dTa' Vn 

(3.42) 

where M denotes 1R+ or iM. depending on Fu. Thus Fu(cr(xa), a~^{xa')) satisfies the ho- 
mogeneous PDE, and, since cr IdwhenXa — > x^, Fu(cr(xB), cr~^(xa')) = K oo{xB,Xa') 
Therefore, Fu(o"(xa), cr"-'^(xa/)) is equivalent to Fu(xa,Xa/) since it satisfies the same 
PDE and boundary conditions (assuming uniqueness). □ 



3.2.2 Neumann Kernels 

The Neumann elementary kernel K^f^ will be defined for U embedded in non-compact 
M. 

Lemma 3.6 Making use of the decomposition introduced in (3.36), the Neumann 
elementary kernel is given by 
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(3.43) 



Proof. Clearly K^'' satisfies the inhomogeneous PDE. For the boundary condition, 



VnaKi^\xa,Xa')\^^ = Vne [-R'oo (Xa, X^/ ) + Fu (x«, X^.)] 

= V„giiCoo(Xi3, X„/) - VnaFu(xs, Xa') 

= 



(3.44) 



where ng is a unit normal to the boundary in the direction of the interior of U. The 
second equality follows because, from the construction, VngX{Ta' , ■)\dv contributes 
opposite signs for paths which contribute to iiCoo and Fu respectively. The third 
equality follows from (3.36) recalling that r,^ ^ when Xa ^ x^. □ 



Note that (3.43) can be written 



K\j ^(x„,x„/) = iiroo(xa,Xa') + Fu(cr(x„),cr ^(x„/)) . 
Lemma 3.7 The Neumann boundary kernel is given by 



(3.45) 



)(x„,x^) := e{x{T^^, z),Xb) exp [-S{x{t^^, z))]vn 
where Vn90(x, x^) = 5(x, x^) and 



(3.46) 



0(x, Xb) := < 




1 



for 
for 



X^Xb 
X)^Xb 



(3.47) 



The ordering is with respect to a foliation induced by a Gauss normal coordinate 
system relative to the boundary. 

Proof. The theorem ensures that K^Q^ is a kernel of the homogeneous PDE. Also, 



= K 



dV 



, because the term vanishes on the boundary since r,^ — > 



as Xa — > Xfi. So the boundary condition is VngKg '{xb',xb) = S{xb',xb). □ 



Corollary 3.2 For Neumann boundary conditions Vng<p{'^B) the solution of the in- 
homogeneous PDE can be written (up to a possible constant) 

fl^(x„) = / iif{f^(x„,x„0/(x„0 dxa> + [ Kf\xa,XB)V^g^{xB) dxB . (3.48) 
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3.2.3 Eigenfunctions and Eigenvalues 



Notice that it is not possible (using the theorem) to determine eigenfunctions satis- 
fying [L — X)^ = with vanishing Dirichlet boundary conditions since this would 
require both f — and (p — 0. However, this tack can be used for eigenfunctions with 
non-vanishing Dirichlet boundary conditions, i.e. f ^ 0. 

Eigenfunctions that satisfy {L — X)^ = which vanish on the boundary can be 
represented in terms of the elementary kernel with the choice oi f = and (p = 0. 
Then (3.33) becomes a Predholm Equation of the second kind; 



The allowed eigenvalues can be determined by estabhshed methods. 



4 Conclusion 



The construction of a path integral that is a solution to a general class of linear second 
order PDEs was presented. Specialization leading to path integral solutions of elliptic, 
parabolic, and hyperbolic PDEs with Dirichlet/Neumann boundary conditions will 
be developed in a subsequent paper. 

This work can be expanded in several directions. One could try to generalize the path 
integral to solve higher order or quasilincar PDEs. It would be potentially profitable 
to check if this construction offers any advantages regarding numerical methods. In a 
physics context, the integrals that have been constructed are based on configuration 
space. These can be recast in terms of integrals based on phase space yielding a (per- 
haps) more fundamental or useful representation. Finally, it is important to extend 
this work to Grassmann variables and fields. 

More significantly, it remains to determine existence and uniqueness of the solutions 
for the general class of PDEs considered here. This is no small task. It is hoped that 
the unifying construction presented here will aid in the effort; and maybe afford some 
insight into PDEs in general. 



Acknowledgment I thank C. DeWitt-Morette for helpful suggestions and discus- 
sions. 




(3.49) 
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A Path Integral Properties 



Recall that // is a complex Borel measure and Q{x, x') and Z{x') are assumed to 
be continuous, bounded, and //-integrable functionals. Many properties of T^q zx 
derived in this appendix are a consequence of analogous properties of Jx> d/i. 

Proposition A.l (Lineeirity) The integral operator fx T^e,zx is a linear operator 
onT{X). 

Proof. 

j {aF^ + hF^) {x) Vq^zx = / ^ Z{x') {adfi + hdv) {x') 

= a / Zix') dfi(x') + b f Z(x') duix') 
Jx' Jx' 

= a / Ffj_{x) Vq^zx + b F^{x) Vq^zx . 

J X J X 



(A.l) 



□ 



Proposition A. 2 (Order of integration) The order of integration over X and X' 
can be interchanged. 

Proof. Prom Definitions 2.1-2.3, it follows readily that 



9(x, x ) dij,{x') 



Ve,zx= / iv,(x) Ve,zx 

J X 

^ [ \[ e{x,x') Ve^zx] dii{x') . (A.2) 

X .J X 



□ 



Proposition A. 3 (Pubini) Suppose that X is the disjoint union of two separable 
Banach spaces, i.e., X = Xi U X2 and Xi n X2 = 0. The dual space X' likewise 
decomposes. IfQ{x,x') = Q{{xi,X2), {x[,X2)) = Qi{xi,x'i) ■ Q2{x2,X2) and Z{x') = 
Z{{x[,X2)) — Zi{x[) ■ Z2{x'^, then 



\ Fi^{x)Ve^zx= / F{xi,X2)Ve^^z2X2 

J X J Xi \_J X2 

= / F{xi,X2)T>Q^^ZiXi 

Xo L'^ Xi 



1^ei,ZiXi 

T^02,Z2X2 



(A.3) 
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Proof. 



f F^{x)Ve,zx= [ Z{x')dfi{x') 

= / , / , Zi{x[)Z2{x'2) dn{x\) dnix'^) 

= / / / / ei(xi,x;)e2(x2,x^) 

Jx[ Jx^ JXi JX2 

= / Q{x,x')V0,^^z2X2'Dei,ZiXidid{x') 

J X' J Xi J X2 

= / / / Q{x,x')dii{x')VQ Z2X2'De^,ZiXi 
JXi Jx2 Jx' 



/ / F{xi,X2)Ve^^z2X2 

J X\ L*' Xo 



(A.4) 



The fourth equahty is a consequence of the Fubini theorem for The order of integra- 
tion over X' and X in the fifth equahty can be interchanged according to Proposition 
A.2. □ 



Note that Vq^ z^Xi and T>q^ z2X2 are possibly different integrators depending on ©1, 
©2 and Zi, Z2. 

Corollary A.l For X an infinite dimensional separable Banach space andM a finite 
dimensional manifold, 

[ f = f f . (A.5) 

Jx Jm Jm Jx 



Proof. Definition 2.1 can be used to define a translation invariant integrator for finite 
dimensional X = M" by 

/ e^p{-^w-„\»)-'^Mu',u)un} i^g^ |-V2 ^ Q^p{-^wMu')} (A.6) 

where VFru is a positive definite invertible quadratic form on R". This integrator can 
be extended to manifolds in the usual way. Then if X2 = K" in Proposition A. 3, 
Ix /r" — /m" Ix and by extension, the corollary follows. □ 



Proposition A.4 (Mean value) Version 1: If there exists an x e X such that 
Q{x, •) = Z{-), then 

F^{x) = / F^{x) Ve^zx . (A.7) 

J X 
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Version 2: Let (x) :— J^^xVqzx and define 0({x), •) := Z(-). Then, 

F^{{x))^ I F^{x)Ve^zx. (A.8) 



Proof. Version 1: 



F^{x)= j e{x,x')dii{x') 
= / Z{x')dn{x') 

= / F^{x)Ve,zx. (A.9) 

t/ X 



Version 2: 



F^{{x)):= [ Q{{x),x')di,{x') 

J X 

= j^^ Z{x')dn{x') 

F^{x)Ve,zx. (A. 10) 



□ 



Proposition A. 5 (Change of variable) Let M : X Y and R : Y' ^ X' be 
diffeomorphisms. Under the change of variable x ^ M{x) = y, 

[ F,{y) V^^y = [ F,{M{x)) Ve,zx (A.ll) 

J 1 J X 

where Q and Z are defined in (2.8) and (2.9). 
Proof. Note that 

(F^oM)(x)= f e{M{x),y')dfx{y') 

= J^^ e{M{x), R-\x')) d^{R-\x')) 
= f e{M{x),R-\x'))di^{x') 

J X 

= J Q{x,x') dv{x') 

= F,(x). (A.12) 
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Hence, 



F,{y)'^e;zy= j Z{y')d^{y') 

= / (ZoR-^){x')di^{R-\x')) 

= I Z(x') diy(x') 
Jx' 

= / F^{x)'De,zx 

F^{M{x))Ve,zx . (A.13) 



□ 



It is important to keep in mind that Vq^x and Vqzx are different integrators. 
However, when M{X) — X and M' is nuclear, the two are related by Restriction 2.1. 



Corollary A.2 If M{X) = X and M' is nuclear and F^{x) e J^r{X), then 

[FMl^e,zy = / F^(M(x))|DetM(',)|2)e,z^ . (A.14) 

Proof. By (2.11), 



- / ^ Z{x') dv{x') 

J X 



= |Deti?'(^-.(,,))|Z(x') di^ix') 
= / F.(x)|DetM(',)|P0^^x 

= / F,(M(x))|DetM(',)|Pe,z^ . (A.15) 

J X 

The fourth line follows from the third by Definitions 2.1-2.3. □ 

Viewing T^qzV ^ a form, this result can be interpreted as 

/ F^= f M*F^ . (A.16) 
Jm{x) Jx 



24 



Loosely speaking, this says that the pull-back of the form in 'local coordinates' is 
{M*F^){x) = |DetM(^.)|F^(M(a;)). Alternatively, the corollary can be interpreted as 
a transformation property of the integrator; 

1^e,z{M{x)) = |DetM(',)| V^^^x . (A.17) 

Similar reasoning yields 

P0,z(lnM(x)) = |Det(lnM)'(i„,)| Pe,z(lnx) . (A.18) 

In particular, for M{x) = x + xq where xq is a fixed element in X or M{x) = e'^^°^"a; 
with e G M+, this characterizes translation and scale invariant integrators on and 
Z+ respectively. 

Proposition A. 6 (Integration by peirts) Suppose that e T^{X), then 

f dyF^{x) Ve^zx = - [ F^{x) dyVe^zx . (A.19) 
If y is any fixed point xq G X, then 

Proof. The first equality follows trivially from Restriction 2.2. 
Now, if X I— > X + Hxq, then Restriction 2.2 and Proposition A. 3 give 



0= / d:,,[F^{x)VQ^zx] 

=11 



SFJx) _ ^ , , S ^ 



a;o(t)o?t. (A.21) 



Then (A.20) follows because Xo{t) is an arbitrary function and Corollary A. 2 implies 
F^{x + hxo) = |DetM(^)|-iF^(x) = F^{x) so that d^,F^{x) =0. □ 



Equation (A.20) can be viewed as an infinitesimal characterization of an invariant 
integrator. For example, ii X — then it characterizes a translation invariant inte- 
grator. Note, however, that it holds more generally; it is true whenever |DetM^'^-j| = 1 
for M : Z^ —>■ Z^. The same argument goes through for z — lax with M{x) — e^^"^^°x 
yielding 

= / ^44^ ^©'^^ = / wAr,1^e,zz . (A.22) 

J\n(X) Sz(t) J\n(X) 5z(t) 
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B Integrators 



This appendix contains a discussion of the well-known Gaussian integrator, the related 
Dirac integrator, and introduces two new integrators; the Hermitian and gamma 

integrators. The Hermitian integrator is a generalization of the Gaussian case, and 
it facilitates integration of polynomials. The gamma integrator is an integrator that 
is inspired by the Laplace transform of a gamma probability distribution in much 
the same way as the Gaussian integrator is inspired by the Fourier transform of a 
Gaussian probability distribution. 



B.l Gaussian and Dirac integrators 
Let 

Q{x, x'; s) — exp {—{Tr/s)Q{x) — 2m {x' , x)-^} 

Z(x'; s) = exp {-7rsW(x')} . (B.l) 

Here s G C_|_, Q is a nondegenerate bilinear form on X = H^{IJ C M) such that 
Rc(Q/s) > 0, and W is its inverse on the dual space X'. Specifically, Q{xi,X2) — 
{Dxi, X2)x and W{x\^x'<^ = {x'l, Gx'2)x such that DG — Idx' and GD — Idx- A one 
parameter family of Gaussian integrators Vq^wx (or T>ujs{x)) is defined according to 
the general scheme by 

j e^^{-i-l^)Q^^)-^^ii^''^)x} Vq^wx := j exp^-''^^<^''^>^> Vujs{x) := exp{-'^^^(^')> . 

(B.2) 

The above characterization is for © : X x X' — > C and Z : X' ^ C For the general 
case e : X X X' and Z : X' ^ C^, (B.2) holds component-wise. 

The Gaussian integrator is normalized, 

VuJsix) = 1 ; (B.3) 
has zero mean, 

/ {x',x)xVuj,{x)^0; (B.4) 

and covariance, 

/ {x[,x)x{x'2,x)xVu;,{x) ^ ^W{x[,x'^) . (B.5) 

JX ZTT 
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In terms of coordinate functions, the covariance can be conveniently written as 

/ x"(t)x'^(u)Pw,(x) = ;fG'"'3(t,u); (B.6) 

where G"'^(t,u) is the (nondegenerate) Green's function of i.e., for x'^ = and 
4 = 5^, W{S^, S^) = (5f , GS^) =: C'^it, u) with a, /3 e {1, . . . , m}. The elements of 
X are pointed paths, and so G'"^(t,u) inherits boundary conditions from the paths. 

The change of variable formula (A. 11) can be used to reduce integrals of the type 
fx f{{x[,x))Vu!s{x) to finite dimensional integrals: 

Proposition B.l Let L:X— ^M" by x ^-^ u = {u\ . . . ,u"} where u* = {x'i,x)x 
for a given x'i e X' and denote W]s_m. = {W o L) so that Qm.^{u) = [W^m]ijVi''u^ : = 
W^l{u), then 

[ f{{x[,x))Vuj,{x)^ [ /(u)|detsiyi,-|-^/2exp{-^^^)^*^'^'"'-'^'<"''">«"> du. (B.7) 

Proof. The proof is straightforward using (A. 11). The determinant factor is a conse- 
quence of the normalization (B.3). □ 
More generally, integrals of the form 

/ f{{x[,x),...,{x'^,x))Vu;,{x) (B.8) 

reduce to finite dimensional integrals and can ultimately be evaluated in terms of G 
as is well known. 

Gaussian integrators possess some important properties that enable the construction 
of path integral solutions to second order partial differential equations. The first is 
the semi-group property ([5]). 

Proposition B.2 Let (t/tF)(x„) := JzdF{xa ■ S{t, z))Vujs{z), then 



{UtM^m^a) = (C/te+t,F)(x„) (B.9) 
where tc > > 0. 

Proof Define the path z : [0, n + rj M'^(C'^) by 

(Ziir) for < T < Th 
- - ' (B.IO) 
ziin) + z2{t - n) for n<T <n^Tc 

where Zi G Z^^ and Z2 G Z'^^- Then 
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= / / F(xa • E{tc, Z2) ■ S{n, zi)) Vu;s{zi)Vu;s{z2) 



(B.ll) 



By the uniqueness of the solution to the differential equation (3.1) and the Fubini 
relation (A. 3), this becomes 

{U-rSUnF)){^a) = / / i^(x„ • S{t, + n, (^1, ^2))) Vus{z^)Vus{z2) 
= / F(xa ■ E{tc + n, z)) Vusiz) 

J 

= (C/U.,F)(x„). (B.12) 

□ 



Next, under an affine map A : X ^ X hy x ^ Ax = Lx + Xq where L is an invertible 
linear map and Xq & X is a fixed element, the image of Vg iyx is 

Vq^wx = \DetL\-^VQ^wA{x) . (B.13) 

This follows from Corollary A. 2. In particular, for L = Id, (B.13) expresses the 
fact that a gaussian integrator is translation invariant. Finally, for F(x) £ J-ji{X), 
gaussian integrators satisfy the integration by parts formula. 

This is just the infinitesimal form of translation invariance of the integrator. 

Related to Gaussian integrators is the Dirac integrator V6{x). The definition of the 
Dirac integrator ([6]) is based on the fact that a delta function can be represented by 
a Gaussian with zero width. It is characterized by 

lim / exp{—{n/s)Q{x) — 27ri{x',x)x}T^Q,wx 

\s\-^OJx 

= : / ex.p{-27ri{x',x)j^}Sx{x)VQ^wx 

J X 

—: / exp{—27ri{x',x)j^}'D5{x) 

J X 

= hm exp {— 7rsVF(x')} 

|s|^0 

= 1 . (B.15) 
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Conversely, the limit |s| — > oo defines the inverse Dirac integrator 

lim / ex.Y>{-{T^/s)Q{x) -2'Ki{x\x)x}'DQ^wx 

\s\^ooJx 

-. f e-x^{-2m{x\x)x}'D5~'^{x) 
= lim exp {—7rsW{x')} 

s|— >oo 

=:6xix'). (B.16) 

In short hand, \im\s\-^oT^^^six) —'■ 'DS{x) and lim|s|_»oo =: 'DS~^{x). 

This definition is 'good' because, under the hnear map L:X— >]R"byxi— >u = 
{u^,...,u"} where u' = {x'i,x)x for a given x'i e X', the integral reduces (using 
Proposition B.l) as expected for any n; 

/ exp{—27ri{x',x)Y}1^S{x)^lim / exp{ — (7rs)lVi^Ju — 27ri(u', u)Rn} 

\detsWRn\-^^^du 
— / exp{— 27ri(u', u)]Rn}(5"(u)du 
= 1 . (B.17) 

Similarly, 

/ exp{—2'Ki{x',x)y-}T)5~^{x)^ lim / exp{— (7rs)H^]^Ju — 27ri(u', u)k"} 



exp{— 27ri(u', u)]Kn}du 
= 5"(u') . (B.18) 

The Dirac integrator possesses the expected properties; 
Proposition B.3 

/ F^(x)V5(x) ^ F^(0) , (B.19) 

and for non-trivial arguments, 

f F,{x)V5{M{x)) = Y: |DetM^^)|-iF^K) . (B.20) 

Proof. First, 
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/ F^{x)V6{x) 

J X 




im / exp{—(27r/s)Q{x) — 27ri{x',x)^}'DQ^wxdiJ,{x') 

I *0 J J X 




(B.21) 



and second, 



/ F^{x)V5{M{x)) 

J X 



[ F^{x)S{M{x))Vq,wx 

J X 

[ F,{M-\ymy)VQ^{M-\y)) 

J Y 

I^F,{M-\y))5{y)\DetM^^]'\VQ^^{y) 
Y,\DetMl,J-'F,{xo) 



(B.22) 



where M{xo) =0. □ 

In short-hand notation, 5{M{x)) — Y.^^ \DetM'^^^-^\~^5{x — Xg). 



B.2 Hermite integrator 

Another integrator of particular interest is a generahzation of the Gaussian integrator. 
For this appUcation, the space of paths is = H^{1] C M). is equipped with a 
quadratic form Q{z) = {Dz,z)zd and the dual space Z^ is equipped with a related 
quadratic form W{z') = {z', Gz')zd' such that DG — Id^d' and GD — Idzd. Assume 
W is nuclear and -Fju(x) e J^r{X). 

Let 



©^(z, z'; s) := Det 




( ^Jt^Q{z)/s\ exp 



ifn(i)exp 



z2_27rj(z',2)} 



{-KQ{z)/s-2TTi{z',z)} 



(B.23) 
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with s e C+ and Hn{z) the n-th order functional Hermite polynomiaL If the space 
has the structure Z — ZiU Z2 ■ ■ ■ U such that Zi n Z2 ■ ■ ■ H Z^ — 0, require 

Hn{zi, Z2,..., Zr, s) = Hn{zi; s)Hn{z2; s) X ... X Hnizr, s) . (B.24) 



Let 



where 



Z^{z'; s) := 1^^, H^- [^7rsW(y'); exp— ^(^'-^') (B.25) 

(B.26) 



Hn'{z'; s) := Hn{z; s] 



Then the n-th order Hermite integrator is characterized by 



/ exp{ ^""'^^ '^^^«} X>p„,s(^) := / Hniz;s)expi 



-P-2-Ki 



Vq^WZ = Zn{z'] S) . 

(B.27) 

It follows from (B.25) and (B.27) that the Hermite integrator is normalized according 
to 

1 forn = 

VpnM = { . (B.28) 



otherwise 



Clearly, VpQ^s{z) is the familiar Gaussian integrator. Integrable functionals are of the 
form 

FJz; n, s) = H^{z) exp"^' / exp{~'^'^^''^^^'} dn{z') . (B.29) 

Under the linear map L : Z'^ ^ Why z ^ {z' z) e M, the quadratic form becomes 

(g o L-'){{z\ z)) =: Q^iiz', z)) = {z', zfjW^z') . (B.30) 
Hence, Proposition B.l yields 



/ Vpn,,{z) 
Jzi 



TISW{Z') 



n/2 



X 



\sW{z')\-^''' 

LHr.{^,u)e.^{-'/^-^i^'ndu. (B.31) 



More generally, let L : Z^ M™ by ;z u and denote W^m = {W o L), then 
(5Km(u) = [l^RijijU^u^ := W^l{u) and 
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fi{z',z))Vp^,,iz) 



det 



ttsWm 



n/2 



I det sWum\-'^^'^ 



X 



(B.32) 

where Hn{u) := 11™ ^ifQ,.(u*) with \a\ :— E^^(ci;j) = n is an m-fold product of Hermite 
polynomials (see e.g. [26]). 

The Hermite integrator facihtates integration of normal ordered (Wick ordered) mono- 
mials: 

Proposition B.4 Define the normal ordered functional monomial -.z"^: by -.z"^: — 

m/2 

Det 



Hm{^TiQ{z)/s). Let L : by z ^ {z', z), then 

[ ■.{z',zr: Vpn,s{z) = \7rsW{z')rn\Snm 



(B.33) 



where : {z\z)™' : is a normal ordered monomial. 
Proof. Under L : ^ R, 



{z',zy 



nsW{z') 



m/2 



Hr, 



( sW{z') '^)) 



(B.34) 



The integral follows immediately from (B.31) and the definition of : z^ :. In a Fock 
representation \n) for the simple harmonic oscillator with ground state |0), 



{z\zy 



nsW{z') 



m/2 



[a + a' ) 



(B.35) 



The binomial expansion and normal ordering yield (see e.g. [27]) 
: {z'.z)"^ : \n) 



TISW{Z') 


rn/2 g / m j 


2 


i=0\i j 



(B.36) 



□ 



Corollary B.l 

^ W : e<-'.^) : Vp^,,{z) = e-^(^') . (B.37) 
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B.3 Gamma integrator 



Specialize to the case where the space of paths Z+ is the Sobolev space i?^(E), i.e. 
z G = if^(E). In order to more clearly differentiate between elements of and 
Z^, it is convenient to now switch notation and write r G = Z^. (Recall that 
T = Inx with X G Z^.) is a Hilbert space with Hermitian inner product 

(Tl,T2)ffi = / Tir2* dt Ti,T2 G T+ . (B.38) 

is 

The range of r is taken to be the path C+ := t(S) = [0, r{tiy)] U [0, T*(tft)] C C+ with 
an associated direction determined by the order of the intervals. 

The dual space T+' is the space of linear forms r' : — > C by 

T ^ (r', t)j,+ := T'{t)T{t)dt . (B.39) 
In particular, define the linear form : — > C by 

T^LV [ f(t)dt , (B.40) 

where a; G C. T+' is equipped with complex Borel measures 

Using the general scheme, characterize a two parameter family of integrators by 

e(T,T';u;) = exp{r^(T) - {T',T)r^+} 



(B.41) 



Z{t'; u, v) = (r' - ; Re(T' - a;) > 

where r G r+ is the path r : E ^ C+, r' G r+', v G C+, and 



_ /^r-fdt _ /.(^)r-rfs 
^ - /,fdt - r(E) ■ ^^-"'^ 



The gamma integrator Di^r (or T>^^^i,{t)) is defined by 



(r' - (vY 
(B.43) 



Note the normalization 



V^.Ar) = i-u^)-" ■ (B.44) 
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The associated space of integrable functionals Tr{T^) contains elements of the form 

F^(r;w) := exp{T^(T)} / ^, exp{-(r', r)^+} ^/.(t') . (B.45) 

On Tr{T^), the integral operator /^+ Vj.t is defined by 

[f^{t-u)V,t := [ J¥-u;)-^dpi{r') . (B.46) 

According to Corollary A. 2 and the characterization (B.43), P^r is scale invariant; 
i.e., in short-hand notation D,y(er) = V„t. This can be seen directly from (B.43), 
which also implies Vj^^^{eT) = Vj^^^^^It). 

Typical applications lead to integrals of the type 

lF^{T{h)-uj)V,T . (B.47) 

where T{ii^ is usually either pure real or pure imaginary. The change of variables 
formula (A. 11) allows this type of integral to be reduced to a one dimensional integral. 
Consider the special case where E = T := [ta, tj C M, —uj e M_|_, and C+ = M+. Then, 
under the hnear map L : T+ — > M+, the change of variable formula (A. 11) and the 
definition (B.43) give, 

/ 2^7a, u{t) = C [ exp{u;u}u^-^ d\x (B.48) 
Jt+ ' Jw+ 

where the constant C is determined by the normalization condition (B.44) to be 
Hence, 

Proposition B.5 Under the linear map L : T+ M+ the gamma integrator reduces 
to 

!v^.Ar)^j ^#^u^#nu). (B.49) 
Jt+ M+ V{y) 

Up to a factor of {—ujY, the integrand of (B.49) is just the gamma probability dis- 
tribution which justifies the appellation 'gamma' integrator. 

Corollary B.2 When T(t) is pure imaginary, the linear map L : T+ iM. induces 
the reduction 



f ^ / ^ 1 n^^ exp|a;u| „ , , 1 f 



=1/ 

2 Jil 



exp{a;u} 

iR r(i/) 



-ioo exp{ci;u*} 

r(^) 



U^'dilliU*) 



u^dilnu) 



(B.50) 



where u 



-u . 
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As another example, let E = R+ and t(E) = with Fn periodic Fn{T{t)) — F^(T(t + 
2n)). 

Proposition B.6 Under the map L : — > S^; the gamma integrator reduces to 

/,I'7„.(r)^(l-e")/"5^!^u-'du (B.51) 

Jt+ Jo -i[v,-uj) 

where 7(z/, —uj) = u'^^—uY iFi^u, u + l;u). 

Proof. The proof follows the same reasoning as Proposition B.l. The periodicity of 
accounts for the different domain of integration and the normalization terms. □ 



Proposition B.7 The gamma integrator satisfies 

riUyV^.Ar) - I, ^7.,(.+.)(r) (B.52) 

for g e C+ and r(S) = C+. If t{T) — then the constant factor on the right-hand 
side is ^{v + —u!)/j{i', —cu). 

Proof. The change of variable r i— > r(tfe) = u G C+ yields 



exp{a;u} 

c+ r(z/) 



u^+^d(lnu) 



r(z/) Jc^T{iy + g) 

r(^) ^ 

+ g) f ^ . . 



(B.53) 



□ 



For g a real integer this gives the moments of the gamma integrator, but note that it 
holds for g G C+ in general. 

Proposition B.8 Under scaling with respect to lo, 

V^e^Ar) = e-" j^^ V^^Ar) ■ (B.54) 
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Proof. The equality follows trivially from the characterization (B.43). □ 

A particular element of the family V^i^^i, that will be useful later is defined by 

[ F{t)Vt:= lim f F{t)V^^,,{t) . (B.55) 

Hence, for integrands of the form F{{t',t)), 

[f{{t',t))Vt^H f F(t)d(lnt) (B.56) 

where C_|_ C C+ is the range of r. The normalization factor M replaces the inde- 
terminate factor exp°°/r(0+) and is fixed by a suitable normalization condition on 
F(t). 

Proposition B.9 Vr is scale invariant and normalized. 

Proof. Scale invariance follows immediately from (B.43) and (B.55). The normaliza- 
tion is 

Ivt^ lim / P7,,,(t)= lim(-^)-^=l. (B.57) 

JT+ v-^0+ JT+ 1,^0+ 

|a;|— >oo >oo 

□ 



C Vciriational Principle 



A key step in constructing path integral solutions to PDEs with finite boundaries 
is to introduce a time reparametrization denoted by r. The first exit time^^ men- 
tioned in the introduction, which is related to r, is well motivated from a stochastic 
point of view. The concept of first exit time and Ito's formula lead readily to the 
functional integrals of interest in the stochastic formulation. Unfortunately, it is not 
immediately obvious how to translate these constructions into a path integral formu- 
lation. Moreover, since the path integral formulation would be expected to include 
the stochastic results as a subset, it would be advantageous to develop the path in- 
tegral constructions based on concepts basic to the path integral formulation rather 
than to simply translate stochastic concepts. Therefore, it is instructive to formulate 

""^^ The term "time" is used here to refer to the evohition parameter of some path. Though 
expedient and conventional, its use can be problematic in some situations; I wiU be more 
careful when warranted. 
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a physical motivation for introducing the parameter r. In doing so, some insight and 
physical intuition is gained for the solutions that are constructed. 



First, recall Feynman's reasoning ([9]) that produced a path integral solution of the 
Schrodinger-type parabolic PDE. Roughly speaking, the prescription is to: choose a 
wave function that represents the state of a physical system at a given point at a given 
time; weight this wave function by the exponentiated classical action functional for 
paths with fixed initial and final positions; and "integrate" over all such paths. From 
this prescription, it is not immediately obvious how to proceed when the physical 
system is contained in a bounded region (See however [10], [11], [12] and references 
therein for an alternative approach.). 

A hint about how to proceed appears when one constructs a path integral repre- 
sentation of a fixed-energy transition amplitude ([6]). The variational principle for 
paths with fixed energy differs from the case of unrestricted paths in that a time 
reparametrization is required ([13]). The fixed energy path integral yields a solution 
of the elliptic Dirichlet problem, and the time reparametrization can be re-interpreted 
as a path reparametrization such that all reparametrized paths which start at a given 
point at an initial time reach some boundary (not some point) at a final time. 

Now, some of the paths of interest for a physical system with a boundary are paths 
that have a fixed initial point and intersect a given boundary at some time. This 
description corresponds to a variational problem from a fixed initial point to a man- 
ifold in the dependent-independent variable space (see e.g. [14]). It turns out that 
this type of variational problem incorporates a variable endpoint in the functional 
integral which can be interpreted as parametrizing the time it takes a path to reach 
some boundary, i.e. a first exit time. Alternatively, the variable endpoint can be inter- 
preted as a non-dynamical dependent variable which takes into account the implicit 
constraint induced by a physical system with a boundary. 

To formulate the variational principle for paths taking their values in a manifold M 
that intersect a boundary, consider the dim(m + 1) dependent-independent variable 
space N = M X M with a terminal manifold of dimension {m+ 1) — k defined by the 
set of equations {5'jk(x, t) = 0} where A; < m, x e M, and t e [ta,tb] C R. Let 



be the functional to be analyzed with a; : M — > M. The extrema of I{x) solve the 
variational problem for point-to-boundary paths. I will refer to the solutions as critical 
paths and denote them by Xcr{t). For the case of M = R", the variational problem is 
solved by the usual Euler equations supplemented by 'transversality' conditions (see 
for example chapter six of [14]). 

There are two limiting cases of interest. When the terminal manifold in N coincides 
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with the boundary in M, then k — 1 and the transversahty conditions reduce to 



F{tb,x{tb),x{tb)) 



uVS{x{tb)) ■ x{tb) 



(C.l) 



where i/ 7^ is a constant. Thus critical paths, which satisfy the Euler equations and 
this constraint, solve the variational problem for paths with fixed initial point that 
intersect a boundary. For free motion, (C.l) implies the critical paths intersect the 
boundary transversally. 

The other case of interest is when the manifold in the dependent-independent space is 
"horizontal", i.e. x(tb) is fixed and the terminal manifold is a line along the t direction. 
The terminal manifold is determined hj k = m equations and the transversality 
conditions yield 



where e is a unit vector in the x direction. If, in particular, F — L + E where L is the 
Lagrangian of an isolated physical system and E is a, constant, then this is just the 
fixed energy constraint {dL/ dx^)x'^ — L = E. Consequently, the variational problem 
in this case is solved by paths with both end-points fixed that have constant energy. 
This gives some insight into the connection between fixed-energy path integrals and 
the elliptic Dirichlet problem. 

To summarize, the presence of a boundary in the variational problem requires the 
introduction of a dependent variable r. For paths ending anywhere on a boundary, 
the critical paths satisfy a transversahty condition. For paths ending at a fixed point 
within a bounded region, the critical paths satisfy a constant energy constraint. Inas- 
much as the critical paths dominate in a path integral, the variational results motivate 
the use of r in a path integral construction. Moreover, the evaluation of a path integral 
often requires explicit expressions for the critical paths: therefore, the transversality 
and fixed energy constraints can be expected to play a role in actual calculations. 
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